Accepted for Publication in Astrophysical Journal Letters 



First Science Results From SOFIA/FORCAST: Super- Resolution 
Imaging of the S140 Cluster at 37/im 

Paul M. Harvey^, Joseph D. Adams^, Terry L. Herter^, George GulP, Justin Schoenwald^, 
Luke D. Keller^, James M. De Buizer'^, William Vacca^, William Reach"^, and E. E. 

Becklin^'S 

ABSTRACT 

We present 37/im imaging of the S140 complex of infrared sources centered 
on IRSl made with the FORCAST camera on SOFIA. These observations are 
the longest wavelength imaging to resolve clearly the three main sources seen at 
shorter wavelengths, IRS 1, 2 and 3, and are nearly at the diffraction limit of the 
2.5-m telescope. We also obtained a small number of images at 11 and 31/im that 
are useful for flux measurement. Our images cover the area of several strong sub- 
mm sources seen in the area - SMM 1, 2, and 3 - that are not coincident with any 
mid-infrared sources and are not visible in our longer wavelength imaging either. 
Our new observations confirm previous estimates of the relative dust optical 
depth and source luminosity for the components in this likely cluster of early 
B stars. We also investigate the use of super-resolution to go beyond the basic 
diffraction limit in imaging on SOFIA and find that the van Cittert algorithm, 
together with the "multi-resolution" technique, provides excellent results. 
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Introduction 



Sharpless 140 is a relatively diffuse H II region at the edge of the much denser L1204 
mol ecular cloud that harbors several clusters of young B stars at a distance of about 900 
pc ( ICrampton fc Fischerlll974j ). The most luminou s object in this cloud was first observed 
i n the near- and mid- infrared by iBlair et al.l ( 119781 ) and by iHarvey. Campbell. &: Hoffmann 
( 1l978l ) in the far-infrared. Blair et al. found a very red, compact object with a spectral 
energy distribution (SED) similar to that of the BN object in Orion in addition to several 
other less remarkable infrared sources. Using the Kuiper Airborne Observatory, Harvey et 
al. found strong, nearly spatially unresolved emission in the far-infrared from 35 - 175/im. 
They fit several crude models to the SED that showed that the dust optical depth to the 
central source had to approach unity in the 50-100/im spectral region and that the total 
luminosity was equivalent to that of an early B star. 

In the intervening years a host of much higher angular resolution near-infrared and radio 
observations have been published. These have revealed many additional infrared sources 



many other molecular clouds harboring young stars, e.g. ( 


Evans, \ 


Vlundv, Kutner & 


DePoy 


1989 




Kraemer et al. 


2001; 


de Wit et al. 


2009 




Preibisch & Smith 


2OO2I 


Kurtz et al. 


1999; 



Bally et al.ll2002l ). It also s eems likely that one or rnore of the strong mid -ir sources contains 



two or more young B stars (jPreibisch fc Smithll2002l : iPreibisch et al.ll2001l) and that the radio 
jet from IRSl is likely to be due to an ionized equatorial wind ( Hoarel 2006 ). Very cold, dense 
condensations only emitting longward of A ~ 300/im ar e also seen and may be even younger 
stage s of star formation within the molecular cloud ( iMinchin. Ward- Thompson &: White 
1995h . 



Most of the luminosity of very young stellar objects and protostars is emitted in the 
far-infrared, wavelengths that are essentially invisible to ground-based telescopes due to ab- 
sorption by water vapor. Until the recent commissioning of the Herschel Space Observatory 
(Pilbratt et al. 2010) and the Stratospheric Observatory for Infrared Astronomy (SOFIA; 
Gehrz et al. 2011), the far-infrared, however, has suffered from a substantial deficit in an- 
gular resolution compared to neighboring wavelengths observable with large ground-based 
telescopes. High mass young stars typically form in dense clusters, and this puts a spe- 
cial premium on the value of high angular resolution to resolve individual components to 
understand their evolutionary state. 

We present here the first imaging observations of the S140 cluster at wavelengths long- 
ward of the ground-based atmospheric window at A ~ 20 - 25/xm with an angular resolution 
of a few arcseconds. With this resolution we are able to distinguish the three main luminos- 
ity sources seen at shorter wavelengths as well as diffuse emission between the objects. We 
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describe the observations in §2] and then in ||3] describe some efforts to improve the angular 
resolution with various deconvolution techniques. In ^we compare our 31 and 37/im maps 
with those at shorter and longer wavelengths and show the spectral energy distributions for 
the individual members of this cluster, and in IIS] we summarize our results. 



2. Observations and Data Reduction 

Our observations were made with the FORCAST camera (Herter et al. 2012) on board 
SOFIA on the night of UT 1 Dec 2010. The FORCAST camera has two separate 256x256 
infrared arrays that cover the wavelength range from ~ 5 - 40/im in a variety of bands 
with 0.768x0.768" pixels. The two arrays can observe simultaneously through a dichroic 
beamsplitter that divides the wavelength range shortward and longward of 26/xm, or the long 
wavelength array can be used by itself without the intervening dichroic. The observations 
presented here were made with the 37.1/im filter (hereafter, 37/im) without the beamsplitter 
and at 11.1 and 31.5/im with the beamsplitter (hereafter "11" and "31"/im). 

In addition to the observations of S140, a number of other science investigations were 
undertaken that are also reported in this issue, and several calibration targets were observed 
for both flux and point-source-response (PSF) calibration. The FORCAST images were 
calibrated to a flux density per pixel using a standard instrument response derived from 
measurements of standard stars and solar system objects over several flights (Herter et 
al. 2012). The color correction from a fiat spectrum source to that of S140 is negligible 
(< 1%). The estimated 3a uncertainty in the calibration due to variations in fiat field, 
water vapor burden, and altitude is approximately ±20%. The point source response has 
been derived from averages over several flights as described in detail by Buizer et al. (2012). 
We discuss below the relatively minor uncertainties involved in using this average PSF, and 
tests of its applicability using fi Cep data taken on the same flight as S140. The nominal 
diffraction limit of SOFIA at 37yum, 1.019A/D with its 14% central obscuration, is 3.1", and 
our observations show that achieved resolution is quite close to this limit, ~ 4" FWHM. 
At shorter wavelengths the SOFIA PSF is affected more by telescope pointing jitter and 
seeing, with s eeing being dominant short ward of 5/im because of the shear layer just above 



the telescope (lErickson fc Dunhamll2000l ). 



Our images of both S140 and fi Cep were obtained with standard thermal infrared 
techniques of chopping and nodding. Because the cluster is relatively compact, we were 
able to chop and nod completely on the ~ 200 arcsec square array. The chop frequency 
used was 2 Hz, and chop amplitude of 90" with a nod amplitude also of 90" perpendicular 
to the chop direction. Each observation resulted in 120 seconds of on-chip integration as 
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the result of one ABBA nod pair. We obtained 10 useful observations of S140 and two 
of /i Cep at 37/im, and two of S140 at 11 and 31/im. The individual observations were 
co-added by cross-correlating all the subsequent images with the first. The sky rotation 
was essentially negligible between the individual S140 observations. The telescope software 
to enable accurate absolute pointing was not yet operational for these observations, so we 
used the highest resolution, nearest w avelength observations (24.5/im) from a large ground- 
based telescope (Ide Wit et al.l 120091 ) to establish the absolute coordinate system for our 
observations. For the purposes of deconvolution discussed later, the S140 frames were left in 
the original instrument coordinate system to match those of the PSF calibration image. 

Figure [1] shows the results of the processing described so far for S140 as well as related 
images from the literature discussed below. Using aperture photometry we have derived 
flux densities for the three obvious components of the image as well as for two additional 
local flux maxima. These results are listed later in Table [1] for a 9" diameter aperture. 
We discuss these further in §1] after we describe our deconvolution work and the resulting 
higher resolution image. The statistical uncertainties in these fluxes are negligible because 
of the very high signal-to-noise ratio, well over 100, for all the sources except for the faintest 
objects at ll/zm where S/N ~ 10. Tabled] also lists the luminosities of these "ob j ects" for 
the assumed distance of 960 pc, using the 50-100/um flux densities of iLester et al.l (Il986l ). 



3. Deconvolution and Super-Resolution 



With NASA's previous airborne observatory, the KAO, iLester et al.l (119861 ) for exam- 
ple showed that significant improvements in the apparent angular resolution in far-infrared 
images were possible through super-resolution techniques such as maximum entropy de- 
convolution. With the KAO's instrumentation these techniques were limited to 1-D scans 
rather than 2-D images. In the intervening years, the improvement in detectors and in 
analysis techniques suggests the possibility of substantial resolution improvements in 2-D 
far-infrared images with SOFIA. To that end we have experime nted with a variety of algo 



rithms and parameters avail able as part of the "MRA" package (jStarck. Pantin fc Murtagh 



2OO2I : IPantin fc StarcklEggsl ) to improve the resolution in our S140 image. Our best results 



were obtained using the van-Cittert algorithm with multi-resolution support for regulariza- 
tion. With this mode the final deconvolved images did not depend strongly on the number 
of iterations, and therefore this was used for all our results below. For example, the final 
image shown used 500 iterations, but the image changed very modestly after as few as 20 
iterations. Additionally, we did not force positivity in the image. 



As mentioned above we have used the PSF derived by Buizer et al. (2012) and done 
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several tests to give us confidence in its applicability to the S140 images. First we deconvolved 
the observed fi Cep image and compared the result with the m uch higher angular resolution, 
but shorter wavelength 24.5/xm image of Ide Wit et al.l (|2008[ ). We found the spatial extent 
and spatial asymmetry to be very similar. We also found the amplitude of the wings to 
the northeast and southwest to be relatively larger at 37/xm, which is entirely consistent 
with the expected dust temperature of order 150K at a radius of a few arcseconds. In a 
second test of our derived PSF, we produced several similar PSF's with small changes in the 
width in both telescope axes, elevation and cross elevation, and then deconvolved the S140 
image with these modified PSF's. In all cases the major details of the deconvolved images 
were essentially the same as for the derived PSF. So we are confident that none of our 
conclusions depends on the small uncertainties in our knowledge of the PSF. Finally, Figure 
[U shows the o bserved and deconv olved images of S140 at 37/im together with the observed 
Subaru image ( Ide Wit et al.ll2009[ ). It is clear that virtually all the structure seen at the 0.6" 
diffraction limit of Subaru at 24.5/im is also visible at 37/im in the deconvolved image, albeit 
at somewhat lower resolution. With the native 37/im resolution of ~ 4", however, much of 
this structure is smoothed to the point of being barely detectable. The final derived PSF is 
equivalent to the theoretical diffraction-limited PSF convolved with a Gaussian of roughly 
1.5" FWHM, and the apparent resolution of the deconvolved S140 image is ~ 2.4". Figure 
|2] shows the derived PSF along with the observed fi Cep image illustrating the larger extent 
of /i Cep at 37/im. 



4. Comparison With Ground-Based Mid-Infrared and Radio Observations 



We have already noted above the excellent agreem ent between our 37u m deconvolved 
image and the much higher resolution 24.5/tm image of Ide Wit et al.l ( l2009l ). In particular, 
in addition to the good agreement in position of the three main peaks, the agreement in 
morphology of the low level, diffuse emission north of IRS 1 and also to the southwest 
suggests that the deconvolution process has not developed spurious features as might have 
been the case with a poor estimate for the PSF or a poor algorithm. There are probably 
still some residual positional problems in our 37/tm image at the level of ~ 1" due to the 
inherent uncertainties in fixing the aboslute positions by registering to the shorter wavelength 
ground-based images, but the agreement between the two images is remarkably satisfying. 
As a corollary it is also clear that there are no new sources present in the 37/tm data, either 
raw or deconvolved, that are not present in the Subaru 24.5/tm image. Most importantly, 
we do not see any evidence for emission above t he background from either of SMM 1 , SMM 
2, or SMM 3, the 450/im sources discovered by iMinchin. Ward-Thompson &: White! (119951 ). 
At 37/im any emission from these sources must be below a level of ~ 30 Jy to be invisible in 
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our image. This confirms their hkely interpretation as very cold, prestellar condensations. 
Four very faint sources in the 37/im image are indicated in Figure [1] that have counterparts 
in the Spitzer IRAC images in the Spitzer archive, presumably low luminosity young stellar 
objects. 

As noted earlier. Tabled] lists flux densities derived from our 37/im image. We list both 
the values measured using the raw image, as well as from the deconvolved image. Since 
the deconvolution algorithm conserves flux, the total flux densities for both methods are 
identical, but we believe the deconvolved image provides a more realistic measurement of 
the flux densities in small beams required for a reasonable comparison of the relative fluxes 
and colors of the individual components in this region. Table [1] also lists the observed 11 
and 31/im flux densities for a 9" aperture. Figure [3] shows the spectral energy distributions 
for the various local maxima that have been identified previously in this region and which 
are clearly detected in our deconvolved map. This shows graphically the fact, as noted by 
previous investigators, that the coldest mid-infrared condensations in this region are the 
diffuse local flux maxima called "NW" and "VLA 4" (labelled in Figure H]), while the hottest 
SED belongs to IRS 1. IRS 2 and 3 exhibit SED's that are intermediate in color temperature 
with IRS 3 being the cooler source. Quantitatively, IRS 1 dominates the shorter wavelength 
emission with 75% of the total emission at 11/xm, but at 37/im it contributes less than 40% 
with the rest arising from the cooler sources, IRS 2, 3, NW, VLA4, and the surrounding 
diffuse emission. 



In their original KAO study of S140 lHarvey. Campbell. Sz Hoffmann! (jl978[ ) showed that 
the dust optical depth was likely close to unity in the far-infrared based on the resolving 
power possible with that instrumentation. With the substantially higher spatial resolution 
available with SOFIA/FORCAST we can re-examine this conclusion. For example, the ratio 
of 37/im to 53/im flux densities shown in Figure |3] for IRSl implies a color temperature of 
order 85K. If we assume that half of the total flux arises in the central 4" core, this implies 
r ~ 0.5 at 53/im, quite consistent with the original KAO study. 

In order to compare quantitatively the diffuse er aission in our 37/im map with that 
seen at much higher angular resolution at 24.5/im by Ide Wit et al.l (120091 ). it is necessary 
to compare the two ima ges at identical lev els of resolution. This technique was described 
as "beam matching" by iLester et al.l (119861 ). Since the native resolution at 37/im is more 
than a factor of 10 worse than that at 24.5/im it is sufficiently accurate to simply convolve 
the 24.5/im image with the 37/im PSF. Figure H] s hows the orig i nal, u ndeconvolved 37/im 
SOFIA image together with the 24.5/im image of Ide Wit et al.l (120091 ) smoothed to what 
would have been observed with our 37/im SOFIA PSF. This shows that there is substantially 
less extended diffuse emission at 24.5/im. Figure|l]also shows the relative 24.5, 31, and 37/im 
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fluxes along a one dimensional cut through IRS 1 from northeast to southwest along the line 
shown in the figure. This comparison shows quantitatively the higher level of extended, 
diffuse emission at the longer wavelengths, consistent with a radial temperature gradient 
away from the center of heating. This efi^ect can also be seen in the more steeply rising 
24.5-37/im SED slopes noted earlier for the sources labeled "NW" and "VLA4" well away 
from IRSl. 

In order to fully model the various components of this compact cluster, it will be neces- 
sary to obtain spatially resolved images at the peaks of their SED's, i.e A ~ 70/im. For the 
3 brightest sources, this will be possible with SOFIA using the upcoming HAWC far-ir cam- 
era. To understand the structure of the cores of the sub-millimeter sources would, however, 
require the sensitivity and higher angular resolution of Herschel with its PACS and SPIRE 
cameras. 



5. Sumrriciry 

We have demonstrated the utility of the FORCAST camera on SOFIA for mapping 
compact regions of star formation at wavelengths beyond the familiar ground-based mid-IR 
atmospheric windows. This allows us to define the SED's and colors of the major emission 
sources in S140 with unprecedented angular resolution in the far-IR. Since the data were 
taken with quite high S/N, we were able to successfully pursue spatial deconvolution of the 
image and show that the basic structure which appears at A ~ 25fim with sub-arcsecond 
resolution from large ground-based telescopes can also be identified in our far-IR images. 
We find emission from all the major sources seen in this cluster at shorter wavelengths from 
ground-based telescopes and confirm that the coolest sources are IRS3 and the diffuse knots, 
VLA4 and NW. We also see no evidence for emission from the sub-mm peaks in this region. 
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Table 1. Aperture^ Flux Densities'^ (Jy) of Sources in S140 



Image 


IRS 1 


IRS 2 


IRS 3 


NW 


VLA 4 


Total^ 


Observed 37yum 


2176 


485 


669 


465 


381 


6730 


Deconvolved 37/im 


2624 


577 


710 


438 


358 


6730 


Observed 11 


110 


4.0 


9.7 


3.5 


2.6 


145 


Observed 31/im 


1585 


368 


401 


242 


184 


3780 


Luminosity (L0) 


10,700 


3200 


2900 


800 


650 


19,150 



'^For 9" diameter aperture. 

^Summed over entire map shown in Figure [H 

'^Absolute calibration uncertainty = 20%; negligible statistical un- 
certainty. 




Fig. 1. — Top row (left to right) - 24.5/im Subaru image of Ide Wit et al.l ( 120091 ). decon- 
volved 37yum image from this paper, and observed 37/im image from this paper; Bottom 
row (left to right) - observed ll/im and 31/im images from this paper, and 450/im images 
from the SCUBA archive (htt p: / / cadcwww.dao.nrc.ca/jcmt / search/product ) with the po- 
sitions of IRS 1,2, and 3 marked with "X", as well as the positions of the sub-millimeter 
peaks SMMl-3. North is up and east to the left. The locations of 4 faint sources that are 
coincident with Spitzer IRAC sources are indicated by arrows in the upper right panel. 



-6-4-2 2 4 6-6-4-2 2 4 6 
Arcsec Arcsec 



Fig. 2. — Left - Derived PSF as described in the text. Contours are at levels of 2.5, 5, 10, 
15,..., 90, and 95% of the peak. Right - Observed image of the nominal PSF calibrator, // 
Cep, with contours at the same level. 
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Fig. 3. — Spectral energy distributions of the various sources detected in the S140 region at 
37fim. For the 37/im flux densities we have used the 9" aperture fluxes hsted in Table [T]from 
the deconvolved map. For the 11 and 31/im fluxes from this study, we used the 9" aperture 
fluxes in Table [1] fron i the observed images. At the sho r ter w avelengths we have used the 9" 
aperture fluxe s from lEvans. Mundy. Kutner &: DePoyI (119891 ) at 1-2.5/im, the 3.45/im flux 
for IRSl frorn Blair et al.l (119781 ). and measured 9" aperture fluxes from the 24.5Atm r nap o f 
de Wit et al. ( 20091 ). At 50 and lOO/im we used the fluxes estimated by Lester et al. ( 1986 ) 



for the three major sources. 
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Fig. 4. — Bottom Panel - 24.5/xm image of Ide Wit et al.l ( 120091 ) convolved with 37/im PSF 
(left), and our observed 37/im image illustrating the more extended, diffuse emission at the 
longer wavelength. The diagonal line through IRS 1 marks the 1-D cut along which the 
relative fluxes are shown in the upper panel, showing quantitatively the higher level of ex- 
tended emission at the longer wavelengths. The approximate positions of all the components 
of this source mentioned in the text are indicated in the lower right panel. Upper Panel - 
relative fluxes along the diagonal line below for the observed 37/im image (dash-dot), the 
observed 31/ini ima ged convolved to the 37/im resolution (dash), and the 24.5/im image of 



de Wit et al.l (120091 ) convolved to the 37/im resolution (solid). 
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